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ABSTRACT 

We present a comparison of the predictions of perturbation theory for the hierarchi- 



cal J-order correlation amplitudes Sj = £ j /£ 2 



J 



10, with the results of 



o 



large numerical simulations of gravitational clustering. We consider two different initial 
power spectra of density fluctuations, one being flatter than the other. The clustering 
amplitudes Sj measured in these models are different at all scales. In each case, the 
perturbation theory predictions give an excellent agreement with the simulations on 
scales for which the variance is approximately linear, i.e. for scales on which £ 2 ~ 1- 
For cells of radius R > 30 /i -1 Mpc, the sampling variance arising from the finite size 
of our simulations dominates the results. We also find that Sj are roughly independent 
of time, £1, or A, as expected. On comparing with the observations, one can use these 
results to discriminate between models for structure formation. 



INTRODUCTION 



The statistical properties of a density field can be charac- 
terised in a simple and efficient way by the J-order connected 
moments (or volume-averaged J-point correlation functions) 

S^J g j ( R ) =< 8 R > c of density fluctuations 8r averaged in 
spherical volumes of radius R (for definitions of the nota- 

Q^tion used in this Letter, see Gaztanaga 1994). 

**^ Gravitational evolution from an initially Gaussian den- 
s ity field generates higher order correlations, { j with J > 2, 
I which are given by the so called 'hierarchical' constants: 

O - " ~ J ~ 1 



Sj 



The values of Sj for density fluctuations 



t/3 



in spherical cells have been estimated using perturbation 



_j theory (PT) in the limit f 2 — > (e.g. Juskiewicz et al. 
1993, Bernardeau 1994b). In this limit, the Sj are inde- 
pendent of time and they are a function only of the shape 
of the initial variance £ 2 . Baugh, Gaztanaga & Efstathiou 
(1994, hereafter BGE94) have performed a detailed system- 
atic comparison of the predictions of PT with the results 
of iV-body simulations, for models with an initial two-point 
correlation, f 2 given by the standard cold dark matter sce- 
nario (SCDM), i.e. for a universe with Q = 1 and h = 0.5, 
specified by T = fi/t = 0.5 (h denotes Hubble's constant in 
units of 100 h kms -1 Mpc -1 ). Using simulations with differ- 
ent numbers of particles and different computational box 
sizes, BGE94 confirm the conclusions of previous studies 
(e.g. Juskiewicz et al. 1993, Kofman & Bernardeau 1994, 
Bernardeau 1994a) and extend them to higher orders. The 
main result of these analyses is that the PT predictions for 
the form of the Sj in the SCDM model agree with the values 
measured from the simulations for scales where f 2 < 1. In 
this regime, according to PT, the values of Sj should be in- 
dependent of time, also in agreement with what is found in 
simulations (BGE94). These results hold for other shapes of 
f 2 and the values of Sj are almost independent of the values 



of the density parameter Q or the cosmological constant A 
(Bouchet et al. 1992, Bernardeau 1994a). 

These results are important because they indicate that 
PT can provide us with predictions for the properties of a 
density field undergoing gravitational evolution which are 
independent of Q, A or time. In contrast the increase in the 
amplitude of the 2-point function, £ 2 , due to gravitational 
evolution does depend upon time, Q, A and the biasing of the 
light distribution with respect to the mass. Hence there are 
several possible factors that could contribute to the 2-point 
correlations, so that a comparison with the clustering in 
galaxy catalogues would not necessarily distinguish between 
them. In particular, it would be difficult to argue that we 
have conclusive evidence that fluctuations grow by gravita- 
tional evolution by just using the observations of £ 2 . Higher 
order moments Sj can help to remove some of this ambigu- 
ity. They also provide us with a powerful tool to study how 
galaxy fluctuations trace matter fluctuations (see Gaztanaga 
& Frieman 1994). 

In this letter, we compare the PT predictions with esti- 
mations made from large iV-body simulations for two models 
which have different correlations at all scales. For com- 
parison, we have also made artificial point distributions with 
an identical size box, number of particles and values of £ 2 , 
and roughly the same number of clusters and mass function 
in clusters. We use this control sample to test if, in realis- 
tic but non gravitational models with clumps, there is also 
correlation between the values of Sj and the shape of £ 2 . 



2 SIMULATIONS 
2.1 N-body models 

The initial linear theory power spectrum of the density fluc- 
tuations used in the simulations is that of Bond & Efstathiou 
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(1984) which is a function of T = fi/t. This form for the 
power spectrum applies for scale-invariant CDM universes 
that have low baryon densities, Qb ~ 0.03. We have simu- 
lated two models that we label T = 0.5 and T = 0.2. The 
T = 0.5 model corresponds to the SCDM model, Q = 1 
and h = 0.5, while the T = 0.2 model is a spatially flat low 
density CDM variant, with a cosmological constant A = 0.8, 
h = 1 and Q = 0.2 at the present epoch. Although these 
are CDM motivated models, what is relevant here is that 
the shapes they predict for the two-point statistic, P(k) or 
f 2 , are different. The results we obtain would also apply for 
other models that produce similar shapes for f 2 . 

We ran ensembles of 5 simulations for each model, using 
the same random phases, in a box of size Lb = 378 ft -1 Mpc 
containing N = 126 3 particles. The simulations are evolved 
using the P 3 M code of Efstathiou et al. 1985. 

The volume averaged correlations { j are estimated from 
moments of counts-in-cells as in BGE94, averaging the re- 
sults over the simulations in each ensemble and using the 
dispersion between members of the ensemble to estimate 
the sampling error. Different stages in the evolution of the 
simulations are labelled in terms of the linear PT variance 
in spheres of radius 8/t -1 Mpc, denoted ag, calculated from 
the initial power spectrum. After a number a of expansion 



our simulations the initial value of ag is a' s 
final value is <Tg 



0.10 and the 

1.0 so that they are evolved over 10 ex- 
pansion factors. Nevertheless, the time resolution for Sj is 
effectively about 3 expansion factors as higher order corre- 
lations are initially affected by a transient arising from the 
the original displacements of the particles from a grid made 
using the Zeldovich approximation (see BGE94). 



2.2 Clumps model 

For comparison we also make a set of "mock" point dis- 
tributions based on a combination of the clustering mod- 
els discussed in Peebles (1980, §61). Clusters are made of 
spherical clumps which are assumed to have N g galaxies (or 
matter points) with a number density varying as a power of 
distance from the clump center. To introduce a hierarchical 
structure, clumps are placed at each step of a Rayleigh-Levy 
random walk. Starting from a random position in a periodic 
box, the next clump in the sequence is placed in a randomly 
chosen direction at a distance drawn from another power law 
distribution. Each clump has a different number of galaxies 
with a "mass function" which roughly imitates that of the 
simulations. 

Of course, this model is not a true representation of 
the physical process underlying galaxy formation. We have 
chosen it because it is convenient and seems to reproduce 
some of the properties of gravitational clustering. In partic- 
ular, one can adjust the parameters in the model to obtain 
a given shape and amplitude for £ 2 . We have run several of 
these models but shall concentrate on ensembles that repro- 
duce closely the T = 0.2 shape for £ 2 . We have generated two 
sets of 5 of these models for each ensemble corresponding to 
different amplitudes of f 2 in a box of size L B = 378 A" 1 Mpc 
containing N = 126 3 particles. Half a slice from the central 
part of one of these simulations is shown in Figure 1, where 
it is compared with slices in the T = 0.2 and T = 0.5 models. 
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Figure 1. Slices of 3 X 189 X 378 Mpc, in the center of the 
simulated box in the clumps, T = 0.2 and T = 0.5 models. 



3 RESULTS 

3.1 Clustering in N-body models 

In linear PT the evolution of the density field is completely 
specified by the 2-point correlation which, for = 1, grows 
with the scale factor a as f 2 (i?, t) = a 2 ^ 2 (R, to), where R is 
the comoving scale and to is some initial time. The linear 
growth is illustrated in Figure 2, which shows the values of 
f 2 for two expansion factors in the simulations. The errors 
are always smaller than the size of the symbols. The solid 
(dashed) lines show the linear PT prediction ^ for the T = 
0.5 (T = 0.2) model, obtained from a numerical integration 
of the initial power spectrum, 



V 
2tt 2 



dkk 2 P(k)W 2 (kR), 



(1) 



where W(kR) is the Fourier transform of the spherical win- 
dow with radius R. Each line is normalized using the ex- 
pansions factors a obtained from the simulation. The most 
noticeable discrepancy due to nonlinearities in Figure 2 ap- 
pears at scales R < 10 h' 1 Mpc where the N-body results 
give larger density fluctuations than the linear PT predic- 
tion. Note however that although the T = 0.2 model has 
less power at small scales in the initial conditions it has 
more power than the T = 0.5 model at all scales for later 
output times. 

The hierarchical scaling predicts = Sj f 2 where 
the Sj are roughly constant. Bernardeau (1994b) has esti- 
mated Sj using PT for spherical fluctuations (top-hat win- 
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Figures 3-5 show the amplitudes Sj in the <r& = 1 
output time of the simulations, estimated from the ratios 
■ in au figures, we also show the predictions from 
PT as solid lines for the T = 0.2 model and dashed lines for 
the T = 0.5 model. Errors corresponding to the dispersion 
over 5 simulations are larger in T = 0.2 model. On com- 
paring the values of { j for individual simulations one finds 
that they do not fluctuate around the mean as a function of 
R but can be significantly shifted at all scales with respect 
to the mean (see BGE94). These overall shifts are produced 
by density fluctuations on scales comparable to the box size 
which are not properly represented in each individual sim- 
ulation. According to this interpretation we would expect 
the T = 0.2 model to have larger sampling errors as it has 
bigger density fluctuations at large scales (Figure 2). 



Figure 2. The variance £ 2 a * two different epochs, erg = 0.40 
and erg = 1, estimated from N-body simulations: models T = 0.2, 
open symbols, and T = 0.5, closed symbols. The solid (dashed) 
lines show the linear PT predictions for the T = 0.2 (T = 0.5) 
model. Dotted lines correspond to clumps models with erg ~ 1.2, 
upper line, and erg ~ 0.5, lower line. 
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Figure 3. The hierarchical amplitudes Sg = C3/C2 111 * ne ^ = 0.5 
CDM model (filled triangles) compared with the T = 0.2 CDM 
model (open triangles). The solid (dashed) lines show the linear 
PT predictions for the T = 0.2 (r = 0.5) model. Open (filled) 
circles correspond to the erg = 0.5 (erg = 0.2) clumps model. 



dow) in terms of the derivatives of the linear f 2 : 

For example, for J = 3 and J = 4, Bernardeau finds: 

o 34 

Z>3 = — + 71 

60712 62 7 2 2 
^ 4 = -1 + -7TH + t7i + T72- (3) 



Thus, S j are not constant as the shape of £ 2 , and therefore 
7 j, changes with scale. For our comparison we have set high 
order derivatives jj = for J > 2 (see BGE94). The above 
values correspond to Q = 1 and A = 0, but the dependence 
on Q and A is negligible (Bouchet et al. 1992, Bernardeau 
1994a). 
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Figure 4. The hierarchical amplitudes Sj = Cj / C2 f° r ^ = 
4 — 6 in the T = 0.5 (filled triangles) and T = 0.2 (open triangles) 
models, at two different output times, erg = 1.00 (left panel) and 
erg = 0.40 (right panel). 



In Figure 4 we also show the <jg = 0.40 output time for 
comparison. Note that the results are very similar in both 
output times, indicating that the values of Sj do not evolve 
significantly with time. 



3.2 Clustering in the Clumps model 

The clumps simulations are constructed to reproduce the 
variance of the T = 0.2 N-body model. Dotted lines in Figure 
2 show f 2 for two ensembles (each with 5 simulations) with 
different amplitudes: <jg ~ 1.2, upper line, and <jg ~ 0.5, 
lower line. The values of S3 are shown in Figure 3 as open 
((78 ~ 1.2) and closed (ag ~ 0.5) circles. These values of S3 
are significantly different to the values in the N-body models. 
Discrepancies are even larger for higher orderss: e.g., at R ~ 
8 ft -1 Mpc the clumps model gives S7 ~ 370 compared with 
S 7 ~ 56000 in the T = 0.2 model. 
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Figure 5. Amplitudes Sj = (j/( 2 for J = 7 - 10 in T = 0.5 
(right panel) and T = 0.2 (left panel) N-body simulations. 



4 DISCUSSION 

For both the T = 0.2 and T = 0.5 models, PT gives excel- 
lent agreement for 5j = 3 i ... i io with the simulations on scales 
where the variance is approximately linear, i.e. on scales 
for which f 2 < 1. More precisely, this scale seems to be 
R ~ 7 A" 1 Mpc in the T = 0.5 model and R ~ 14 A" 1 Mpc 
in the T = 0.2 model. For scales R > 30 ft -1 Mpc, the sam- 
pling variance arising from the finite size of our simulation 
boxes dominates our measurements of the Sj. The accuracy 
of the agreement does not seem to deteriorate too much 
when the order increases up to J = 10. 

In Figures 2 and 4 we have compared the results at two 
different output times, separated by 2.5 expansion factors. 
The results indicate that the values of Sj do not evolve much 
with time. This is in spite of the fact that { j evolves by a 
factor ~ 10 J_1 between as = 0.4 and as = 1. 

Visually, in Figure 1, the T = 0.2 model looks closer to 
the T = 0.5 than to the clumps model, but nevertheless f 2 
for the clumps model is closer to f 2 for the T = 0.2 model, 
by construction. This illustrates the importance of higher or- 
der statistics. It is not difficult to construct artificial point 
distributions which have similar 2-point statistics to the sim- 
ulations but very different higher order correlations. For in- 
stance, one could place all particles inside a thin wall within 
the simulated box to produce a given f 2 (see Gaztanaga & 
Yokoyama 1993). Here we have tried more realistic attempts 
in order to see if the PT predictions are a general feature of 
models with "clumps" or a specific feature of gravitational 
clustering. Figure 3 shows how different the values of S3 
for the T = 0.2 gravitational model are compared with the 
skewness in a clumps models with a similar f 2 . We have also 
generated clumps models with different shapes and ampli- 
tudes of f 2 but have not been able to reproduce the values of 
S3 in either the T = 0.2 or the T = 0.5 models. This is not 
surprising, as a visual inspection of the clustering pattern 
in the maps (c/ Figure 1) shows strong differences between 
gravitational and clumps models: clusters are not spherical 
in the N-body simulations. 



The values of Sj provide a way to discriminate between 
different models of structure formation: realistic distribu- 
tions of points can have significantly (within the sampling 
errors) different values of Sj at all scales. In particular, we 
have illustrated here how two distributions with similar val- 
ues of the 2-point statistics f 2 can produce very different 
values of Sj. The distributions that are evolved under grav- 
ity seem to match the PT predictions. In particular, there 
is a correlation between the shape of f 2 and the values of 
Sj, which we do not find in the mock clumps models. More- 
over, these predictions are practically independent of Q and 
A (Bouchet et al. 1992, Bernardeau 1994a), in agreement 
with what we found in the N-body simulations. 

The values of Sj in the APM galaxy distribution 
also seem to follow the PT predictions (Gaztanaga 1994, 
Gaztanaga & Frieman 1994). To interpret this results one 
has to consider first the problem of biasing, i.e. how galaxies 
trace matter. After considering the posibility of both non- 
linear and non-local biasing, Gaztanaga & Frieman (1994) 
have found that the simplest interpretation for the values of 
Sj measured in the APM is that there is little or no biasing 
on scales R > 7 ft -1 Mpc, as otherwise one has to invoke 
a delicate, model-dependent balancing of terms so that the 
estimated values of Sj follow, at each order, the PT predic- 
tions. This interpretation implies that the observed values 
of Sj in the APM are compatible with the clustering that 
emerges from gravitational growth of small (initially Gaus- 
sian) fluctuations, regardless of the cosmological model we 
assume for the universe, i.e. Q, A, Ho or the nature of dark 
matter. 
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